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ABSTRACT 

We present the highly significant detection of a quasi-periodic flux modulation with a period of 22.2 min 
seen in the X-ray data of the Sgr A* flare of 2004 August 3 1 . This flaring event, which lasted a total of about 
three hours, was detected simultaneously by EPIC on XMM-Newton and the NICMOS near-infrared camera 
on the HST Given the inherent difficulty in, and the lack of readily available methods for quantifying the 
probability of a periodic signal detected over only several cycles in a data set where red noise can be important, 
we developed a general method for quantifying the likelihood that such a modulation is indeed intrinsic to the 
source and does not arise from background fluctuations. We here describe this Monte Carlo based method, and 
discuss the results obtained by its application to a oihtr XMM-Newton data sets. Under the simplest hypothesis 
that we witnessed a transient event that evolved, peaked and decayed near the marginally stable orbit of the 
supermassive black hole, this result implies that for a mass of 3.5 x 10^ M©, the central object must have an 
angular momentum corresponding to a spin parameter of a ^ 0.22. 

Subject headings: black hole physics — Galaxy: center — Galaxy: nucleus — accretion — X-rays: observa- 
tions — Methods: data analysis 



L INTRODUCTION 

Detecting Keplerian motion in the accretion flow orbiting 
around a black hole (BH), in the case where an event, local- 
ized to a portion of the flow, retains its coherence over sev- 
eral orbits, entails being able to detect periods of the order 
of lO^'^s (or frequencies of -lOkHz) for solar mass objects. 
The problem is entirely different when considering a super- 
massive black hole (SMBH), and in particular Sgr A*. For a 
non-rotating BH of 3.5 x 10^ M©, the period at the marginally 
stable orbit (MSO, r(Mso) = Ir^ - 6GM/c^) is 1592 s or about 
26.5 min. So as much as detecting kHz quasi-periodic oscil- 
lations may be a delicate matter for certain reasons, as we ex- 
plore long periodic signals, low frequency noise and window- 
ing effects become important, and must be taken into account 
in a proper characterization of the signal. 

In general, to test for a periodic signal, we must know what 
is the probability distribution function (PDF) of the test statis- 
tic in the absence of such a periodic signal, so that a reli- 
able significance can be assigned to the claimed detection. In 
Fourier analysis, a PDF of the form e^^, where Z is the Fourier 
power gives the classical probability measure. For time series 
with important fluctuations or flares (colored noise), one usu- 
ally obtains a distribution that is broader than e"^ resulting 
in overestimated significances. Similarly, windowing effects 
introduced by the testing of non-integral periods with respect 
to the total observation time also affect the probability distri- 
bution. These effects can be taken into account in order to 
derive an accurate probability for a given peak in the power 
spectrum. 

For Sgr A*, the MSO — where we expect the minimum 
fundamental Keplerian period — oranges from Pmi,, ~ 300 to 
^max ~ 1600 s, for a mass of 3.5 x 10^ Mq and spin parame- 
ter > a = T^ < 0.999. Here, J is the angular momentum, c 
is the speed of light, and G the gravitational constant. A num- 
ber of X-ray observations of Sgr A* have revealed flares with 
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durations between a 3 000 and 10000s. We adopt a standard 
minimum quality factor, Q, of 4-5 cycles to claim a period- 
icity, and thus find that for the longest period of 1600 s, the 
minimum flare duration must be 8000 s. The number of In- 
dependent Fourier Spacings (IFS or trials) in our full interval 
(8752 s) is given by: IFS = T,^,,(l/P^i„ - 1/P„,,) = 24. 

Ascenbach et al. (2004) claimed the detection of 5 peri- 
ods ranging from 100 to 2250 s each of which were identified 
with one of the gravitational cycUc modes associated with ac- 
cretion disks. The work of these authors was based on X-ray 
data from XMM-Newton. Genzel et al. (2003) claimed the 
detection of an ^ 17 min period detected during two near-IR 
flares from Sgr A*. Clearly, a confirmed detection of a peri- 
odic or quasi-period signal present in a flare from the central 
BH would be very useful for constraining the nature of the 
flares, and more importantly, as evidence for the geometry of 
the space-time around Sgr A* through indirect BH spin mea- 
surements. In the context of this paper, we use the term quasi- 
period to refer to the average period of a signal that exhibits a 
periodic quality, while possibly comprising an evolution from 
longer to shorter periods, which would leads to a broader peak 
due to a spread in the power over several frequencies. 

2. OBSERVATIONS AND METHODS 
2.1. Observations 

In the last few years, observations of the central BH 
with the Chandra and XMM-Newton X-ray telescopes 
have revealed s everal flares with di f ferent durations and 
spectral indices jBasanoff et al. 200H iBag anoff et al. 2003t 
iGoldwurm et al. 2003t iPorauet et al. 2003). Recently, two 
more flares were detected by XMM-Newton and their prop- 
erties are discussed by Belanger et al. (2005). 

Here, we restrict our analysis to XMM-Newton data, and 
pay particular attention to two data sets: 0516-0111350301 
and 0866-0202670701. These observations of Sgr A* were 
respectively performed on 2002 October 3 and 2004 August 
31, and each comprises one major flare, that lasted ~3 and 
9 ks respectively. Details on the observation and characteris- 
tics of the first flare are presented and discussed by Porquet 
et al. (2003), and those of the other event by Belanger et al. 



(2005). In addition, we analysed the data of pointing 0789- 
0202670601, which also comprised a flare from Sgr A*, but 
solar proton contamination caused a large number of data gaps 
during the flaring events thus making these data unreliable for 
period analysis. As described in Belanger et al. (2005), the 
light curve can be reconstructed under the condition that the 
binning is large enough to be insensitive to the presence of 
very short good time intervals that would otherwise lead to an 
overestimate of the count rate. In this case, 500 s was found to 
be a reasonable bin size, but this does not allow a meaningful 
period analysis for an event that lasted ~3 ks. 

For each data set, we used the PN event list constructed by 
selecting good X-ray events in the 2-10 keV energy range, 
within a radial distance of 10" from the radio position of 
Sgr A* (Yusef- Zadeh et al. 199 ?). The PN instrument has 
a larger collecting surface, better sensitivity and most im- 
portantly, a much higher time resolution than the MOS 
cameras — 73 ms compared to 2.6 s in FullWindow mode^. 
This makes the PN camera the XMM-Newton instrument of 
choice for timing analysis. Figure^show the light curve of 
the 2004 August 31 flare and of the 2002 October 3 observa- 
tions, both with 120 s resolution. 




LSxIO* 



Fig. 1. — Top panel: Light cm-ve of the 2004 August 31 flare (r~9ks). Bottom 
panel: Light curve of complete - 1 5 ks 2002 October 3 observation. Resolution is 1 20 s. 



2.2. Test-statistic and Periodogram 

The Z^ -periodogram is constructed as follows. For a given 
test period Pj, the phase 0„ of each arrival time f,-, is calcu- 
lated as (pj = intiJP j. These phases are then used to calculate 
the Z^ value defined by Buccheri et al. 1983 as 



T m 



k=\ 



^ cos(kcf>i) + ^ sm(k4>i) 
1=1 ' V/=l 



(1) 



Here, A^ is the number of events and m is the number of 
harmonics over which the sum is performed. We used the 

PN's time resolution can be as high as 0.3 ms in the standard Timing mode with a 
duty cycle close to 100%, and 7^s in the special Burst mode. The duty cycle, however, 
is only 3% in the latter mode of operation. 



Rayleigh statistic (m-l) and thus performed the sum over 
the first harmonic only, since this is the most powerful test for 
sinusoidal signals (Leahy, Eisner & Weisskopf 1983). 

The Rayleigh statistic, Z^ is distributed as a x^ with two 
degrees of freedom (dof). This is so because cos0 and sin0 
are Gaussian distributed, the square of a Gaussian variable 
is c/i/j distributed, and the sum of ;^'^ variables is also a chi^ 
variable for which the number of dof is given by the sum of 
the dof from the individual variables. Thus, the Z,^, statistic 
is x\m distributed, de Jager (1994) showed that the predicted 
significance for a sinusoid from the Rayleigh power is 



■logio(Prob) = 0A34(Np'-/4 + 1), 



(2) 



where A^ is the number of events, and p is the pulsed fraction. 
As mentioned above, the range of periods related to 
the MSO we consider physically interesting is 300-1 600s. 
Nonetheless, periodograms are shown from 100 s to Tob^/3. 

2.3. Background 

In the simplest case, the background is composed of white 
or Poisson noise: fluctuations for which the range of values 
is consistent with counting statistics. A more general case 
considers the background as a combination of white and red 
noise. Red noise (sometimes referred to as 1// noise) arises 
from the inherent source variability and has a power spectral 
distribution that goes as / ", where / is the frequency and a is 
the slope of the power spectrum in log-log space — the larger 
the amplitude of the variations, the steeper the slope. 

Clearly, it is essential to determine the relative importance 
of the component of red noise with respect to the Poisson 
noise in order to accurately estimate the probability associated 
with a given peak in the periodogram. This can be done in dif- 
ferent ways, one of which was recently suggested by Vaughan 
(2005) and entails determining the index a by performing a 
linear fit on the log-log periodogram and then calculating in 
an analytical way the confidence levels, conveniently repre- 
sented as lines parallel to, and lying above the fitted power- 
law. This method is simple, does not require Monte Carlo 
(MC) simulations, and is practical for data that has an im- 
portant component of red noise and whose duration is long 
enough to allow for a suitably large range of frequencies with 
a well defined power spectrum over which to perform the fit.^ 

Another way to characterize the component of red noise in 
an observation is to construct the power spectrum of the ap- 
propriately smoothed light curve. The power spectrum will 
then be dominated by the features associated with the inherent 
source variability, and will exhibit the features of the underly- 
ing red noise component, which can then be characterized as 
described above. This procedure should give results closely 
compatible with the method of de-trending — fitting the light 
curve with a polynomial and subtracting it — which leaves the 
flattened light curve free of the inherent source variability. 
Here, however, the stochastic nature of the red noise is not 
accounted for. 

We have used a method that allows us to work directly with 
the event list and to take into account all noise components 
present in the data in a natural and transparent manner that 
also takes into account the random nature of the noise. This 
was done with the use of MC simulations to generate event 
lists having the same statistical properties as the data. 

* A proper fitting technique must be used in order to determine the power- 
law index (see Goldstein, Moiris & Yen 2004). 



2.4. Pseudo-Random Event Lists 

Most analyses aimed at estimating the significance of a pe- 
riodic signal use MC simulations. We have used such a tech- 
nique to estimate the probability of the null hypothesis, that 
is, the probability that a peak in the periodogram is caused by 
a background fluctuation. However, several different methods 
can be employed, and thus we detail the one we have used. 

In counting photons from a given astrophysical source or 
radioactive decays from an unstable isotope, each event oc- 
curs independently of the previous but with a certain regu- 
larity given by the average count rate. These processes can 
be described by the same statistical law, characterised by an 
exponential PDF that defines the distribution of time inter- 
vals between two consecutive events and for which the mean 
corresponds to the average count rate. Therefore, in order to 
construct a true pseudo-random light curve of total duration T 
and average count rate r, one must draw T xr numbers from 
an exponential PDF with a mean of r. Each of these cor- 
responds to the Af between two consecutive events. Arrival 
times are calculated by summing these Af, and the light curve 
is finally constructed from an event list that has the same sta- 
tistical properties (Poissonian for r - const), as the data. 

In the case where the mean is not constant, we can gener- 
ate event lists with equivalent statistical properties as the data 
set by drawing the Afs from different exponential PDFs with 
means determined by the sliding average count rate calculated 
from the data over an appropriately sized window. 

The size of the smoothing window determines the range of 
periods accessible, (period longer than this size are smoothed 
out). By performing the simulations using different window 
sizes determined on the basis of the length of the observa- 
tion as r/4, r/8, r/16, r/32, etc. down to a specified min- 
imum, the probability associated with each trial period is de- 
termined taking into account fluctuations in the base level of 
the light curve that occur on timescales of the smoothing win- 
dow size. Therefore, in each range of accessible trial periods, 
the significance of each peak is properly assessed. A period 
will be picked out only if it is significant compared to general 
trends in the light curve. Drawing numbers from PDFs with 
means that follow the variation in the count rate takes into 
account the overall noise properties of the data by simultane- 
ously characterizing the inherent variations in the base level 
of the light curve, and allowing fluctuations about this level 
that are proportional to the amplitude of the variation. This 
is the method we used to estimate the probability associated 
with each value in the Rayleigh periodogram. 

In practice, for a given window size, we determine a set of 
mean rates and errors. For a mean rate r + a-y, a Af is con- 
structed in two steps: one pseudo-random number, r, is drawn 
from a Gaussian N(r,o"r), and another, Af, is drawn from an 
Exponential E(r). The phases are constructed from the ar- 
rival times as described above. The full range of testable peri- 
ods (chosen as lOO-r/3) is subdivided in 100 logarithmically 
spaced period bins, and logarithmically sampled by 150 trial 
periods, Pj, to ensure that each bin contains the same number 
of trials. For each Pj, we obtain a Z- value histogrammed in 
the corresponding period bin. Finally, the probability is de- 
termined directly from the Z^ values by finding the quantile 
inverse for the Z^ value of P j in the data set. 

3. RESULTS 

Figures 12 and 13 show the Z^-periodograms, and their as- 
sociated probability as a function of period derived from the 



MC simulations for the flare data sets of 2004 August 3 1 and 
and 2002 October 3 respectively. The In Fig.0] we present 
the periodogram and probability of the data obtained on 2002 
October 3 before the flare. The vertical scale is the same on 
all three Figures. 
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Fig. 2. — Flare of 2004 August 31: Z^-periodogram and associated probabilities 
derived from MC simulations of 10^ event lists. Total flare duration was ~9000s and 
thus the period range is 100-3000 s. The peak is at /*= 1330 s and rises to a Z^ value of 
46, and has an associated probability of ~10"^. 
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Fig. 3. — Flare of 2002 October 3: Z^-periodogram and associated probabilities 
derived from MC simulations of 500 event lists. Total flare duration is ~3000 s and 
period range spans 100-1000 s. No significant peaks are detected 




Fig. 4. — Data of 2002 October 3 before flare: Z" -periodogram and associated proba- 
bilities derived from MC simulations of 500 event lists. Total duration of observation is 
~ 10000 s and period range spans 100-3000 s. No significant peaks are detected 

As expected, the probability function is perfectly anticor- 
related with the value of Z^ over most of the range. Devia- 
tions are apparent towards to end of the range, where Pj be- 
comes commensurable with Tobs- This effect is most notici- 
ble in Fig. |4] where the last three peaks in the periodogram 
grow in height whereas their probabilities hover around one, 
indicating indicates that the MC simulations indeed model the 
overall statistical properties of the data accurately and thus re- 
produce features associated with the noise. Notice that simi- 
lar behaviour, although not as marked, is also seen in the 2004 
August 31 as well as 2003 October 3 events. 



In Fig. 12] the strongest peak at 1330 s has an associated 
probability of 3 x 10"^. The other two noticable peaks at 
370 s and 850 s have respective probabilitites of 0.006 and 
0.007, and are therefore not significant. Peaks with proba- 
bilities greater than 10"^ cannot be considered significant if 
we apply a 3cr requirement to claim a detection and thus, no 
peaks in the periodograms other than the one at 1330 s can be 
taken as significant. The after-trial significance for this signal 
is 7 = 4.14, and using Eq.|2we find a pulsed fraction of 0.13 
(A^ = 2023): a very reasonable number, 

Fig.|5]presents three phasograms resulting from folding the 
flare light curve with a phase of 1330 s. These were con- 
structed in the 2^ keV, 4-10 keV and 2-10 keV bands. There 
is a clear asymmetry in the overall light curve, and a marked 
peak present at low energies, making the 2^keV band radi- 
cally different from the much smoother 4-lOkeV light curve. 
This peak is shifted by ~266 s from the centre of the sinu- 
soid, corresponding to a distance of ~0.5 AU for light. This 
intriguing feature will be investigated elsewhere. 




Fig. 5 . — Light curve of 2004 August 3 1 flare folded with a phase of 1 330 s, in three 
energy bands: 2-10 keV (top), 2^keV (middle) and 4-10 keV (bottom). 



Finally, we also analyzed the pre-flare period of the 2002 
October 3 observation to search for the periods reported by 
Aschenbach et al. (2004) in these data. The periodogram is 
shown in Fig.|4]and is remarkably flat over most of the range 
of physically interesting periods. As we can see from the as- 
sociated probabilities, there are no significant peaks present. 

4. SUMMARY AND CONCLUSION 

We detected the presence of a 22.2 min quasi-period in the 
X-ray data of the flaring event in SgrA* that occured on 
2004 August 3 1 while simultaneous XMM-Newton and HST 
observations were under way (see Yusef-Zadeh et al. 2006, 
Belanger et al. 2005 and Belanger et al. 2006 for details on 
the results of the campaing). The probability that such a 



signal arises from statistical fluctuations was determined to 
be 3 X 10-6 from MC simulations and is conservatively esti- 
mated to be less than ~7.2 x 10"^ if we consider the 24 IFS in 
the period range 300-1600s ( 1). It is important to clarify that 
although an evolution of the period as a function of time may 
be present (as expected from the spiraling motion), and that 
it can be qualitatively seen in the light curve as alluded to by 
Liu et al. (2006), our statistics do not permit us to claim the 
detection of such a feature. 

We considered and analysed two other XMM-Newton ob- 
servations of Sgr A* using the same method for the probability 
estimation: the one that contains the 2004 March 3 1 flare, and 
the one of the 2002 October 3 flaring event. The former was 
found to be unreliable due to solar proton contamination, and 
no significant signals were detected in the latter, in contrast 
with the results obtained by Aschenbach et al. (2004) using 
the same data set. 

The method we have developed to detect a periodic sig- 
nal, and evaluate the associated probability is well suited for 
searches involving long periods with respect to the total ob- 
servation time, and in which fluctuations in the average count 
rate do not follow Poisson statistics. This MC based method 
is just as powerful in the absence of red noise but clearly not 
necessary since the probability can be computed analytically. 
This straightforward and reliable method for quantitatively es- 
timating the significance of a given periodic or semi-periodic 
modulation will be useful in future applications of this kind, 
and hopefully help, at least in part, relieve some of the ambi- 
guities that arise in this type of analysis. 

The radius of the MSO in units of r^ = GMjc^ is given by: 



[3 + Z2 - V(3-zi)(3+Zi+2z2)], 



(3) 



where z, = 1 -i- (1 - a^yi^[{l + aY'^ + {I - ayi\ and zi = 
(3a^ + Z[)''^. The Keplerian frequency is: 



InGM 



(r^/2 + a)- 



(4) 



Taking a mass of 3.5 x lO^'M© and assuming that the peri- 
odic signal originates from orbital motion at the MSO, the 
1330 s period implies that the BH is spinning at a rate given 
by fl a; 0.22 (in the prograde case, the only considered). 

Thus, the main result presented here can be interpreted as 
the signature of a transient event that evolved over about 1 ks, 
simultaneously giving rise to an X-ray and near-infrared 
flare. The emitting region spiraled around the mildly spin- 
ning SMBH, near or at the last stable orbit. In this simple 
scenario, we would have directly witnessed the accretion of 
matter through the accretion disk in the Kerr metric around 
the central dark mass (see Tagger & Melia 2006). 
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